Pairing Gaps, Pseudogaps, and Phase Diagrams for Cuprate Superconductors 
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We use a symmetry-constrained variational procedure to construct a generalization of BCS to include Cooper 
pairs with non-zero momentum and angular momentum. The resulting gap equations are solved at zero and 
finite temperature, and the doping-dependent solutions are used to construct gap and phase diagrams. We find 
a pseudogap terminating at a critical doping that may be interpreted in terms of both competing order and 
preformed pairs. The strong similarity between observation and predicted gap and phase structure suggests that 
this approach may provide a unified description of the complex structure observed for cuprate superconductors. 

PACS numbers: 



o 
o 

Of 
=5 



-a 
c 

o 
o 



> 

00 
ON 
^1" 
(N 
O 

o 

-i— > 

c3 



i 

C 

o 
o 



X 



I. INTRODUCTION 

The mechanism responsible for high-temperature super- 
conductivity remains unresolved despite intense study over 
the past two decades. The observation of pseudogaps and a 
relatively universal phase diagram are thought to represent 
key aspects of the solution, but no theory describes all ob- 
servations in a simple, unified way. Previously we introduced 
an SU(4) model of competing antiferromagnetism (AF) and 
d-wave superconductivity (SC), used coherent states to pro- 
vide a connection to Landau-Ginzburg and symmetry con- 
strained Hartree-Fock-Bogoliubov theory, and demonstrated 
that SU(4) symmetry implies no double occupancy i 1 ' 2 ' 3 - 4 This 
has defined a new, unconventional method to simplify the 
rather complex many-body problem in cuprates which is oth- 
erwise very difficult to solve. 

The SU(4) coherent state solution represents a generaliza- 
tion of the BCS pairing theory having Cooper pairs with non- 
zero momentum and angular momentum. Dynamics of these 
interacting Cooper pairs at different hole-doping and temper- 
ature yields rich phase diagrams that may explain the main 
features of observations. The occurrence of a pseudogap, a 
critical value of doping, and splitting of the SC gap at low 
doping region are all natural consequences of the dynamics. 

In this article we extend our previous discussion and com- 
pare gap and phase diagrams computed using the formal- 
ism developed in Ref. 4 with available data. We demonstrate 
that our analytically-solvable gap equations with doping- 
independent interactions already describe the basic features 
of experimental gap and phase diagrams for the cuprates. 
With addition of simple doping dependence for interaction 
strengths in the low-doping region, our results can reproduce 
quantitatively the superconducting gap, the pseudogap, and 
the doping variation of both the superconducting transition 
temperature T c and the pseudogap transition temperature T* 
extracted from a range of cuprate data. 



II. GENERALIZED SU(4) COHERENT STATES 

In the SU(4) model the configuration space is built from co- 
herent pairs formed from two electrons (or holes) centered on 



adjacent lattice sites: spin-singlet D and spin-triplet % pairs 1 . 
The 15 generators of SU(4) then consist of two operators for 
singlet pairs (D and D^), six operators for triplet pairs (n and 
7T'), three staggered magnetization operators £2, three spin op- 
erators S, and a number operator h. An SU(4) Hamiltonian 
restricted to one and two-body interactions can be determined 
uniquely, 

H = en- GqD^D -Gi^-n- %M.M+ kS- S, (1) 

where Go, G\, and % are effective strengths of li-wave singlet 
pairing, triplet pairing, and AF correlations, respectively. We 
shall assume the total spin to be zero and ignore the last term 
for the present discussion. 

The SU(4) coherent state 2 - 5 is | ) = ST | 0*), where |0*) 
is the physical vacuum and the unitary operator 2T is a 4- 
dimensional matrix parameterized in terms of variational pa- 
rameters u and v, with u± + v 2 = 1 . It implements a quasi- 
particle transformation on the D-% pair space that preserves 
SU(4) symmetry (implying no double occupancy 3 ). The cor- 
responding matrix elements for one- and two-body operators 
are discussed in Ref. 2 and matrix elements for the operators 
appearing in Eq. ([]]) are easily evaluated in the Q. — > °° limit 4 . 



III. TEMPERATURE-DEPENDENT GAP EQUATIONS 



Introducing the "gaps" 



representing correlation energies for singlet pairing, triplet 
pairing, and AF, respectively, the variational Hamiltonian 
H' = H — An (with chemical potential A) is 

(H') = (e -X)n- (A^/Go + A 2 /^ +A 2 /*) . 

Variation of (H 1 ) with respect to u± and v± yields 

2u±v±(e± - A) - A±(u\ - v 2 ) = 0, 

which has the solution 



(2) 
(3) 



4 = £ii 



E± - A 



e± 



e± 

[( £± -A) 2 - 



4 = iH- 



e± — A 



(4) 



-A ± 2 ] 



21* 



: A,, ± A„ 



2 



where we define 

£± = £TA,. (5) 

From Eq. (0]), the definitions for Ad, A K , and A q , and a theoret- 
ical doping fraction x = 1 — n/£l with n the electron number, 
we obtain the temperature-dependent gap equations 

G £2 



which exists only for x less than a critical doping 



A,/ = —. — (w + A + + w_A_ ) 
(w+A + — w_ A_) 
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(6) 
(7) 



4A 

^ = w + (A, + A')+vv_(A, / -A') (8) 

-2x = w+{A q +X')-w-{Aq-X') (9) 
where A' = A — £ and 

l-n±(T) 



w± 



e± 



(10) 



The quasiparticle number densities n±(T) are assumed to be 4 

2 

(11) 



n± 



l+ezp(Re±/k B Ty 

where n± = and w± = l/e± if T = 0. The gaps and A fol- 
low from the algebraic equations ©-Q, the energy E can be 
calculated from (fJJ and 

E={H') + nX, (12) 

and e±, u±, and v± are determined by Eqs. (0]). This defines a 
complete formalism permitting calculation of general observ- 
ables. 

These results are similar to the BCS theory but with impor- 
tant differences. Here we have two pairing energy gaps and 
two kinds of quasiparticles, implying complex behavior (ulti- 
mately tracing to the non-abelian commutator algebra for the 
SU(4) operators 3 ). The quantities e± correspond to two sets 
of single-particle energies {e±}, split by 2A q in the AF back- 
ground. Each level can be occupied by one electron of spin 
up or down. The corresponding pairing gaps are |A±| and the 
probabilities for single-particle levels to be unoccupied or un- 
occupied are u±, and v±, respectively. When G\ =% = K = 0, 
Eq. ([T} reduces to a pairing Hamiltonian, A K = A q = 0, and 
Eqs. (O-dUi reduce to the usual BCS equations. 



IV. SOLUTION OF THE GAP EQUATIONS 

We now describe the solution of the gap equations ©-(O, 
first for temperature T = and then for finite T. Data sug- 
gest that in hole-doped cuprates the interactions in ©I-© 
are attractive, with AF correlation strongest and triplet pair- 
ing weakest. Assuming that % > Gq > G\ > 0, we find the 
general T — solution 



XO.[(x- l -x)(x g -xp 2 



A d = \G Q.[x{. 



^-*)] 1/2 

1/2 



A K = ^G D.[x(x q -x)] 
A' = -\xQx q {l-Xqx)-\G\£lx, 



(13) 
(14) 
(15) 
(16) 



X-Gi 



1/2 



(17) 



The dependence of x q on three elementary interaction 
strengths shows explicitly that the critical doping point results 
from competition between correlations. A trivial A q = A K = 
(pure singlet) solution exists also in the entire physical doping 
range < x < 1 , 



Ao = A d = {G Q a[{l-x z )} 



2N11/2 



A' = -lG Qx. (18) 



Trivial solutions valid for < x < 1 that correspond to pure 
triplet pairing, pure AF states, and metallic states (all A, = 0) 
may be found also, but insertion of A q , A c [ and into Eq. © 
indicates that the nontrivial solution is always the ground state 
where it exists (0 < x < x q ), and that for x > x q Eq. ( fT~8b gives 
the ground state. Thus, when x > x q the ground state is a pure 
singlet pairing state given by ( fT8l , but when x < x q the solu- 
tion (IT3T>— (fT6b (which differs from the singlet pairing state in 
permitting finite values for A d , A K , and A q ) is the ground state 
and ( TT8l defines an excited state. The critical doping point 
Xq separates these two qualitatively distinct ground states and 
marks a quantum phase transition. 



V. CUPRATE GAP DIAGRAMS 

We now use these results to construct gap diagrams for 
cuprate superconductors. The parameters Go, G\, and % are 
effective interaction strengths within a truncated space. In Fig. 
QJ, we first make the simplest possible approximation and as- 
sume that they are constants, independent of doping. Four 
energy scales are shown: 

1. Aq, which measures AF correlations; 

2. the singlet pairing gap A c / [Eq. (fT4bl, which is the su- 
perconducting gap for x < x q ; 

3. the singlet pairing gap Ao = A<j [Eq. (TT8ll, which is the 
superconducting gap for x > x q but is not the ground- 
state order parameter for x <x q ; 

4. the triplet pairing gap A K . 

The AF gap A„ is maximal at x = 0, decreases rapidly with 
doping, crosses the pairing gaps, and vanishes at the critical 
doping Xq. Like A q , the triplet gap A K exists only within the 
doping range < x < x q ; it peaks at x q /2. For x > x q the 
singlet gap is the curve Ao but below x q it splits into two curves 
(A ( [ and Ao) having different doping behavior. Also shown in 
Fig. QJ are gap data from Refs. 6 - 7 that support the complex 
gap structure suggested by the SU(4) quasiparticle solutions. 
In particular, the predicted splitting of the singlet pairing gap 
and termination of the A q gap at a critical doping P q ~ \x q are 
consistent with the data points plotted. 

FigureQJ suggests qualitative agreement between predicted 
and observed gap structure at T = but there are two quan- 
titative discrepancies: (1) the slope of A q is too small at low 
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FIG. 1: (Color online) Energy gap diagrams at T = (a, c) and 
phase diagrams (b, d) for hole-doped cuprates. The doping rate P 
is defined as P = (£1 — n) j£l e with Sl e and £2 being the number 
of lattice sites and the maximum allowed number of holes, respec- 
tively. In (a) and (b) constant interaction strengths are assumed: 
G\ = 3.74, Go = 8.2, and % = 13, while in (c) and (d) doping- 
dependent strengths indicated by the insets to the figures are used. In 
all plots P q = 0. 18, except for the dotted line in (d) marked T q , which 
illustrates how the phase boundary marked T q shifts if P q = 0.16. 
The green dashed curve in (b) and (d) is the empirical T c curve: 
T c = 7c. max [l — S2.6(P — 0.16) 2 ]. Data in (a) and (c) are taken from 
Refs. 6,7 , and in (b) and (d) from Refsi 8 ] 9 ] 10 ] 11 ! 12 ] 13 , respectively. Note 
that data in (b) and (d) marked with the legend (Dai) taken from Ref. 8 
and include data taken from Refsi 11,12 ' 13 . 



doping, and (2) the onset of a pairing gap at zero doping in the 
calculation contradicts data suggesting that this onset should 
occur closer to 5% hole doping. These discrepancies indicate 
that the constant coupling strength assumption gives unrealis- 
tic strong pairing and too weak AF strength in the low-doping 
region. They can be removed by allowing a simple doping 
dependence for the coupling strengths according to the fol- 
lowing physical considerations. First, the SU(4) symmetry 
requires that the pure AF state can exist only when Go = 0, 
which suggests that pairing strength must have an onset at 
P > 0. Second, to reproduce the experimental Neel temper- 
ature, X has t° be strongest at P — 0. These suggest a dop- 
ing dependence for Go, G\, and % given by the inset to Fig. 



QJ, which implies an onset of pairing strength at P = 0.05 (in 
the simplest picture we assume the same onset for singlet and 
triplet pairing, but with different magnitudes), and a large AF 
correlation at half-filling that decreases as doping increases, 
but finally is stabilized as pairing is established. With that 
change, the more quantitative agreement with data shown in 
Fig.QJ results. 



VI. SCALING BEHAVIOR 

In Figs.QJ andQJ: all energies are scaled by £(r c ) max . The 
solution (fT4l indicates that A c [ depends only on x q , once Go 
is fixed. If the critical doping x q has a universal value (for 
example, see Ref. 6 ), the doping dependence for Aj also should 
be universal if gaps are scaled by the maximum T c for each 
compound, since the maximum T c is proportional to Go (see 
Ref£). It is well known that cuprate SC pairing gaps exhibit 
such scaling. 



VII. COMPETING ORDER AND PREFORMED PAIRS 

The emergence of the critical doping x q and the splitting of 
the singlet gap for x <x q are a direct consequence of compet- 
ing pairing and AF correlations below x q . When doping is low, 
AF correlations dominate pairing and a state with large AF 
gap and suppressed pairing is the ground state. Therefore, the 
superconducting gap for the ground state is small (A c / < Ao) 
and the large pairing gap Ao is associated with an excited state 
when x < x q . However, as doping increases the SU(4) sym- 
metry implies that pairing correlations decrease less quickly 
than the AF correlations (see Ref. 1 ) and they eventually dom- 
inate. The energy is then minimized by larger pairing and 
diminished AF correlation. The critical point x q is the hole- 
doping fraction where the ground-state AF correlations vanish 
completely in the variational solution. 

The competing-order picture 6 assumes that the pseudogap 
(PG) is an energy scale for order competing with supercon- 
ductivity that vanishes at a critical doping point where the 
competition is completely suppressed. From Figs, la and lc, 
A q (an order parameter of the AF phase) has precisely these 
properties. But the AF operators are generators of SU(4), 
so A q also is the stabilization energy for a mixture of pre- 
formed singlet and triplet SU(4) pairs that condense into a 
strong superconducting state only after AF and triplet pairing 
fluctuations are suppressed by hole doping. This represents 
a non-abelian generalization of the standard phase-fluctuation 
model 15 for preformed pairs. Thus, the SU(4) pseudogap state 
results from competing AF and SC order expressed in a basis 
of singlet and triplet fermion pairs, which may be viewed as 
a unification of the competing order and preformed pair pic- 
tures. 
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VIII. THE ROLE OF TRIPLET PAIRS 

Triplet pairs are an essential component of the SU(4) many- 
body wavefunction (for example, the SU(4) algebra, which 
imposes the no double occupancy condition on the lattice, 
does not close in the absence of triplet pair operators). How- 
ever, Fig. QJ indicates that the triplet pair correlation energy 
is small in the underdoped region and zero for doping larger 
than P q . The primary role of triplet pairs in the hole-doped 
cuprates appears to lie in fluctuations mediating the AF-SC 
competition. 



unsettled. Our theoretical T* is clearly defined as the temper- 
ature at which the PG in the SU(4) model goes to zero. Once 
the interaction strengths are fitted to the gap data, there are no 
free parameters to modify the theoretical T* curve. Thus, any 
reproducible data that lie substantially above our predicted T* 
curve may be associated either with physics that goes beyond 
the minimal model that is described here, or with a definition 
for T* that is inconsistent with ours. 



IX. CUPRATE PHASE DIAGRAMS 

In Figs.QJ andQJl we show phase diagrams resulting from 
finite-temperature calculations. Figure [TJ) assumes the cou- 
pling strengths of Fig. QJi and Fig. [TJl assumes the coupling 
strengths of Fig. [TJ, so the only adjustable parameter in either 
case is R = 0.6. The empirical factor R < 1 [see Eq. ( fTTT il is 
introduced because actual single-particle energies due to ther- 
mal excitation are non-degenerate and realistic quasiparticle 
excitation should be easier than in our degenerate approxima- 
tion. There are four distinct phases in Figs.QJ andQJl: 

1. An antiferromagnetic phase (labeled AF). 

2. A superconducting phase (labeled SC). 

3. A transitional phase with all three correlations present 
(labeled AF+SC). 

4. A metallic phase. 

The correlations involved in each phase are indicated in paren- 
theses in Fig. QJ> and the doping-dependent temperatures T c , 
T*, and T q (or 7p mark the phase boundaries. Data from 
Refs i 8 i 9 i 10 are compared with the predicted SU(4) phase struc- 
ture in Figs. QJ> and [TJl. The constant coupling approxima- 
tion (Fig. QJ>) reproduces data qualitatively, but the predicted 
superconductivity extends too low in doping and the upper 
boundary on the AF phase does not rise steeply enough at low 
doping (leading to a Neel temperature a factor of 2-3 too low). 
These discrepancies are removed in Fig. [TJl, which allows the 
evolution of coupling with doping indicated in the inset. 

It is commonly believed that in the cuprate phase diagram 
there is a boundary (for example, the one labeled T N in Refii) 
distinguishing an AF long-range ordered state from a disor- 
dered one. We do not discuss this possibility here but note 
that the potential to describe such a boundary is contained in 
the present model. If (contrary to our minimal assumptions 
here) we permit onset of the triplet pairing strength at some- 
what lower doping than for the singlet pairing, a new bound- 
ary will appear in the very-low doping region (P w 0.02-0.05) 
that separates a pure AF phase from the PG region. 

There are additional (generally higher-temperature) data 
that have been related to the pseudogap discussion beyond that 
displayed in Fig. 1. Their physical interpretation is generally 



X. DISCUSSION OF SUPPORTING DATA 

The results of Fig. QJ are supported by various additional 
data. We mention a few representative examples. NMR 
measurements 16 using strong magnetic fields to suppress SC 
in Bi2Sr2- x La l CuOg + 5 concluded that the PG coexists with 
SC and terminates in a critical point near P ~ 0.21. Figure [JJ 
could accommodate P q ~ 0.21 with small parameter changes, 
but also we note that magnetic fields could increase P q by al- 
tering the SC-AF competition. From specific heat, NMR, and 
transport data, Tallon et a)¥. conclude that high-r c phase be- 
havior is universal, with a critical doping P ~ 0. 19 where the 
PG vanishes and where large changes in quantities like the su- 
perfluid density indicate crossover from weak to strong super- 
conductivity, and that this transition implies vanishing short- 
range magnetic order. A key role for magnetic correlations in 
PG states also is suggested by the effect of impurities on c-axis 
optical conductivity 18 . In-plane resistivity measurements 19 
indicate that the phase diagram is surprisingly universal, with 
a PG terminating near optimal doping and different electronic 
states either side of the termination. These measurements also 
indicate a crossover above T c between two kinds of PG be- 
havior (see also Ref. 20 ), which may be explicable in terms of 
the finite but decreasing pair fluctuations expected in the PG 
region, as we now discuss. 

Figure [JJ may explain the vortex-like Nernst signal 21 ob- 
served above T c . The singlet pair gap vanishes there but the 
many-body SU(4) wavefunction has finite pair content be- 
tween the curves T* and T c that decreases with increasing T 
and decreasing doping. Thus, contours for pair fluctuations 
above T c will be similar to observed contours for Nernst signal 
strength, but a Meissner effect is expected only below T c . Re- 
cently, Ong et al 22 concluded that a consistent explanation of 
pseudogap and Nernst data requires PG and SC pairing states 
that are distinct but related by symmetry, as proposed here. 

More detailed properties of Fig. QJmay be investigated in 
future work. For example, the transition between the AF+SC 
and SC phases is 2nd-order if P q lies at higher doping than 
the point b (solid curve T q in Fig. QJl), but is lst-order if it 
lies at lower doping (dotted curve T' starting from the point 
a). This implies a wealth of testable consequences for the gap 
and phase structure depending both on global properties and 
on microscopic details. 
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XI. SUMMARY AND CONCLUSIONS 

In summary, finite-temperature SU(4) coherent states have 
been used to construct a rich, universal gap and phase struc- 
ture for hole-doped cuprates that may be expressed in a BCS- 
like formalism with two kinds of quasiparticles and competing 
superconducting and antiferromagnetic order. We find a pseu- 
dogap of antiferromagnetic character terminating at a critical 
doping P ~ 0. 1 8 that is distinct from the superconducting state 
but related to it by a non-abelian symmetry. The correspond- 
ing pseudogap states may be interpreted in terms of either SC- 



AF competition or preformed SU(4) pairs that condense into 
a singlet <f-wave superconductor as hole doping suppresses 
fluctuations. Our results represent a minimal variational so- 
lution of competing antiferromagnetism and d-wave super- 
conductivity on a fermionic lattice with no double occupancy. 
Therefore, we believe that the general gap and phase structure 
presented here will be a necessary consequence of any real- 
istic theory that takes a doped Mott insulator with competing 
d— wave pairing and antiferromagnetism as the basis for de- 
scribing high-temperature superconductivity. 
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